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Abstract 

We consider the influence of tlie electron-electron interaction on the non- 
linearity of the current- voltage characteristic of the tunnel junction at low 
bias (diffusive anomaly) in the presence of the classical magnetic field. We 
present the theory of a new phenomenon which manifests itself in the strong 
anisotropy of a diffusive anomaly for different orientations of the magnetic 
field with respect to the interface of the tunnel junction. The nonlinear dif- 
ferential tunneling conductance has a universal magnetic field dependence, 
so that only the magnetic field component perpendicular to the interface is 
involved. In particular, when the magnetic field is parallel to the interface, 

the I-V characteristic does not depend on the value of the magnetic field. 
PACS numbers: 73.40. Gk, 71.10.Pm 
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I. INTRODUCTION. 

As it was shown by Altshuler and AronovS, the electron-electron interaction strongly 
influences various physical properties of electron systems. Among those are conductivity, 
spin susceptibility and the one-particle density of states at the Fermi level. Suppression of 
the one-particle density of states at the Fermi level leads to appearance of small deviations 
from Ohm's law in the current-voltage characteristics of a tunnel junction at small volt- 
ages V (zero-bias diffusive anomalies). These authors used the diagrammatic perturbative 
method. A more complete description of the nonlinearities in the I-V characteristics of tun- 
nel junctions which allows to go beyond the perturbative treatment was given in Refs.( 
(Earlier the non-perturbative treatment of renormalization of one-particle density of states 
due to electron-electron correlations was given in Ref. P). It was shown that the strength of 
the diffusive anomaly can be characterized by the effective frequency dependent impedance 
{u ~ eV/h) which the electron feels when tunneling. Immediately after the electron tunnels, 
the system acquires an extra energy due to the interaction between this electron and the 
electrons of the liquid. Therefore, the electron density perturbation must spread in order to 
reach the final state. In this process the many-body system acquires some action which is 
proportional to the impedance mentioned above. The value of this effective impedance can 
be estimated as the resistance of the metal piece which the electron covers in time ~ h/eV 
during its diffusion motion after the tunneling. The value of this impedance governs the 
magnitude of the deviations from Ohm's law. When the conductivity of the system is quite 
large, the effective impedance is much smaller than the quantum unit of resistance h/e^ and, 
as a result, there is only a small deviation from the Ohmic behavior. In the present work 
we will assume that this condition is satisfied. The classical magnetic field influences the 
charge spreading (through the conductivity), hence, the zero-bias anomaly depends on the 
magnetic field value. 

We present here the theory of a new phenomenon which manifests itself in the strong 
anisotropy of a diffusive anomaly for different orientations of the classical magnetic field 
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with respect to the current direction. For definiteness, we considered the tunnel contact of 
two 3D semi-infinite samples and the magnetic field vector can have an arbitrary angle with 
the interface plane (see Fig. |1|). 

As is common, we consider the case of a metallic regime [epTp ^ 1, Tp is the momentum 
relaxation time). We have solved the equations of classical electrodynamics which describe 
the propagation of the charges and electric fields for this geometry and obtained the following 
results. The nonlinear differential tunneling conductance has a universal magnetic field 
dependence, so that only the magnetic field component perpendicular to the interface is 
involved. This is related to the fact that the charge spreading after the tunneling has 
a completely different character for the classical magnetic field oriented perpendicular or 
parallel to the interface. Namely, the magnetic field perpendicular to the interface strongly 
reduces the transverse conductivity a^x of the system (in the region of classically strong 
magnetic fields), hence, the spreading of the charge along the interface plane is strongly 
blocked. The effective impedance increases and has a universal dependence on the magnetic 
field [1 + (cUcTp)^] (multiplied by the impedance value for a zero magnetic field). This means 
that the nonlinear tunneling conductance also has this universal magnetic field dependence. 
For the magnetic field parallel to the interface the charge spreading along the interface plane 
is not blocked because this spreading easily occurs through the Hall drift along the interface. 
As a result, the I-V characteristic for this orientation of the magnetic field does not depend 
on the magnetic field value. This phenomenon resembles the absence of magnetoresistance in 
the metal at low temperatures for one group of carriers. In the case of arbitrary orientation 
only the perpendicular component of the magnetic field is involved because the contribution 
of the parallel component is canceled by the Hall electric fields. 

The size effects in one-particle density of states in the presence of electron-electron cor- 
relations were studied in Ref. 0. This consideration was restricted to the case of a zero 
magnetic field. In contrast, the infiuence of the classical magnetic field on the electron den- 
sity of states was studied only for the case of bulk metali. The authors of the latter did not 
take into account the important role of the surface of a metal in this phenomenon and could 
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not obtain the strong dependence of the zero-bias anomaly on the orientation of an external 
magnetic field. 

The results of present work explain the observations of the recent experiment^. In the 
latter the magnitude of the diffusive anomaly for the electron tunneling between two heavily 
doped GaAs samples was measured for two different orientations of the external magnetic 
field parallel and perpendicular to the interface. It was observed that the magnitude of the 
diffusive anomaly does not depend on the magnetic field value in the case when the magnetic 
field is parallel to the interface and strongly depends on it in the case of perpendicular 
orientation. 

II. FORMALISM. 

Let us consider the concrete geometry of the experiment. We assume that the leads 
are two 3D semi-infinite conducting samples isolated by a thin high barrier of width d (see 
Fig. 0). If a negative voltage is applied to the left lead, the electron may tunnel through 
the barrier from the left lead to the right one leaving a hole behind. The magnetic field 
applied governs the charge spreading process after tunneling. In general, the magnetic field 
is arbitrary oriented so that the diffusion coefficients -Dfcm depend not only on the magnetic 
field value, but also on the angle 9 between the magnetic field and the plane of the insulating 
layer. 

We use the tunnel Hamiltonian method. The total Hamiltonian of the problem consists 
of the following parts: 

H = Hr + Hl + Himp + Ht . (1) 

Here Hr^l describes the electron system in the right and the left leads with the electron- 
electron interaction taken into account. The Hamiltonian Himp describes the electron scat- 
tering by the impurities and the tunnel term Ht = A + describes the transitions of 
electrons through the insulating layer. 
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A= J dh£r'T{r,r')xjj^{r)ij{r') . (2) 

x>d/2 
x'<-d/2 

The corresponding expression for the current operator is 

f = i[A-A^) . (3) 

As usual, the current is calculated in the second-order perturbation theory in the tun- 
nel Hamiltonian. For arbitrary voltage V the current can be obtained0 by the analytical 
continuation of the following Matsubara Green's function determined at uJn = 27iTn , 

ire^'^"" (fA(T)AUO)) 

-eV+iO 



I{V) = 2elm IjdTe''^"^ (ti(r)iT(o) 



(4) 



From now on we assume that e > and ^ = 1. This Green's function is calculated in 
the standard fashion by averaging over the electrodynamic field configuration^ and can be 
expressed in terms of the effective action which describes the low-frequency electrodynamic 
modes of the system excited during the tunneling. When the action is small {S{t) <^ 1), 
one can get the following simple expression for the small correction to the current-voltage 
characteristic at T = 0: 

where Rt is the tunnel resistance of the junction. 

At low frequencies {ujTp ^ 1, is the momentum relaxation time), when the diffusion 
approximation is valid, one can get the following formula for the action!: 

1 ^.3. 



5H = -y rf^rp_^(r)0^(r), (6) 

where 0aj(i") is the electrodynamical potential excited by the tunneling electron, and Pi^(r) = 
/ dr exp{iujT)p{r,T) is the Fourier transformation of the probability density to find the 

— oo 

electron (hole) at point r if Vq is the point of tunneling on the interface of the tunnel 
junction. This quantity is just a common electron-hole (diffusion) propagator which stands 
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in usual diagrammatic approachfl. P(^(r) obeys the diffusion equation in imaginary time in 
the presence of the magnetic field (different signs correspond to different leads), 

l'^IPa.(r) - E Dkm'^k'^mPujir) = sign{Lj)6{r - To) , (7) 

k,m=x,y,z 

with the boundary conditions that there is no diffusive current across the surfaces of the 
leads: 

J2^'-^k'^kPu,iT^) x=±d/2=0 (8) 
k 

This boundary condition is imposed at the distance from the surface which is larger than 
mean free path. It should be noted that in the case of uj^Tp ^ 1 the diffusion coefficients are 
renormalized due to the existance of skipping orbits at the distance of the order of cyclotron 
radius Rc from the surface. However, in the diffusion approximation we consider here the 
contribution of the sliding current to the action ^(ci;) is small with respect to the parameter 
qRc <^ 1, where is the characteristic length scale in the problem.0 

In the coordinate system when the magnetic field is parallel to the 2;-axis, the diffusion 
coefficients have the following form {uc is the cyclotron frequency): 

Dzz = Do = , -Dxx = Dyy = ^ ^ ^2^2 ' (9) 

C p 

Dxy = -Dyx = -sign{uj) Dq ^""^^ 2 ' (-'-^) 

For an arbitrary orientation of the magnetic field in the x-z plane the diffusion coefficients 
can be obtained from the diffusion matrix D (Eqs. (p|,p!OD) by the rotation over the angle 6 
around the y-axis: D{e) = tj^{e)DU{9). 

We assume that the conductivity a of the leads is large, so for the frequencies in question 
{uj ^ o") we use the electroneutrality condition (Debye radius = 0) when calculating the 
electro dynamical potential 0(^(r). In this case there are no charges in the 3D leads except the 
2D charges ±.Qu){j\\) on their surfaces, and (pujij) obeys the system of the following equations 
in imaginary time: 
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J2VkUuJ,r) = 0, (11) 

k 

jk{uj, r) = - ^ o-fcmVm0a;(r) , (12) 

m 

^^'^^^^\\'>\^=±d/2 " - esign{uj)6{rii) . (13) 

Conductivity is related to the diffusion tensor through the Einstein relation, akm = 
e^Vp , where is the density of states for noninteracting electrons at the Fermi sur- 
face. Outside the leads 0(^(r) obeys the Laplace equation. Since for the small frequencies 
considered here the characteristic spatial scale of the variation of 0^; inside the leads is much 
larger than d, we can neglect the electric fields in the leads while finding the relation between 
values of (p^j and at the surfaces. Hence, we obtain a simple relation 

,,=±,/2 = ±27rdQ^(r||)/6., (14) 

where is the dielectric constant of the insulating layer. 

It should be noted that according to Refs.( 0-^) the value Z{lu) = e~'^ujS{uj) can be 
thought as the effective frequency dependent impedance which the electron feels when tun- 
neling. 

Finally, we would like to mention the connection between the method used in the present 
paper and the diagrammatic approachB. Calculation of the correlator in Eq.(^ corresponds 
(in the lowest order) to the account of diagrams a) and b) in Fig. ^.i Diagram a) is the 
exchange correction to one-particle density of states in each lead. Diagram b) contributes to 
the correlator due to excitation of the gapless interface plasmon during electron tunneling. 
We do not take into account Hartree diagram c) which is small0 in parameter Xp/ro <^ 1 
because we consider the case of the weakly non-ideal electron gas. 



III. CALCULATION OF ACTION S{lo). 



Differential Eq. (0) with boundary condition (|^) can be easily solved by the Fourier 
transformation with respect to the y and z coordinates along the surface plane. Then, for 
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the charge density in the right lead {x > d/2) we obtain 

e signioj) 



exp 



iq.D-^D,, + A (x - d/2) 



(15) 



where A'^ = D 



-1 



\UJ\ 



^yy^l + (-^22 ~ ^xx^lz) ■ A similar procedure with Eqs. ([TT| - |l^ 



for the electrodynamic potential 0^ gives 

e sign{u) 



[X, q) 



{2'nd) ^ei\uj\ + da^x - iqyCr.. 



exp 



xy 



{}(lz(^xx^xz + 5){x- d/2) 



(16) 



with 5^ = a^^ 



<^yyTy + \(^z 



a-^a"^ )q^ 

XX XZ / ^2 



. Then, for the characteristic frequency scale of 
the problem {uj ~ -Dg^, q <^ l^^ and U = VpTp) the term (27r(i)^^ei|co'| is small compared 
to others in the denominator of the right hand side of Eq. (|l^), provided d ^ r'^/k (r^ = 
Do/ (To)- The last inequality is usually satisfied for tunnel junctions, and we omit this term 
in the expression for (p^^ in further calculations. This approximation actually means that the 
spreading of the charges in the left and right leads occurs absolutely independently. It is 
determined by the diffusive modes and is not effected by surface plasmon modes which are 
not excited in the frequency range considered (it means that we neglect the contribution of 
diagram b) in Fig. ^). 

By inversion of the coordinates and the sign of the electron charge in Eqs. ( ITsj [16|) we 
can obtain the potential and the charge density in the left lead. Then, S{uj) in Eq. (^) is 
given by the integral over the right lead multiplied by 2. Carrying out the integration over 
X and over q, we present >S'(ti;) in the following universal form: 

^2 



S{uj) = So{u)F{H,d), Soioo) 



2'Kao\/\(jj\D^ 



F{H, 



l + ujlr^sin^e. (17) 



IV. NONLINEAR DIFFERENTIAL CONDUCTANCE. 

In the expression for S{lj) (see Eq. (|13)) the dependence on u is extracted in factor 
So{uj), i.e. the action for H = 0. Factor F{H,6) depends only on the magnetic field. 
Hence, to obtain the correction to the current-voltage characteristic we can easily perform 
the analytical continuation in Eq. (p^ , 
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It seems to be more convenient to represent this result in terms of the correction 6R{H, 9) 
to the Ohmic resistance of the tunnel junction. To do this, we integrate the above equation 
over V , 



so that the correction to the resistance at zero bias is 6R{H, 6) = SR{0)F(H, 6), and 5R{0) 
can be estimated by comparison of the two terms in Eq. (p!9D at eV ~ l/r^ : 6R{0)/Rt ~ 



Function F{H, 9) governs the dependence of the zero-bias anomaly on the magnetic field. 



to the orientation of the magnetic field. In particular, if the magnetic field is exactly parallel 
to the interface, then F{H, 0) = 1 and the correction to the resistance of the tunnel junction 
does not depend on the magnetic field value. In our opinion, this quite interesting effect 
has the same physical origin as the absence of the magnetoresistance in the metal for one 
group of carriers. Namely, additional spreading of the charge along the interface plane due 
to the Hall electric fields that appear in this geometry cancels exactly the blocking of the 
charge spreading due to decreasing a^x in the magnetic field. In contrast, for the case when 
the magnetic field is perpendicular to the interface plane, the Hall electric fields do not 
occur and the charge spreading is governed only by (J^x. As a result, for = 7r/2 we have 
F{H, n/2) = 1 + uj'^Tp, which is just the magnetic field dependence of a^^. It is worth 
mentioning that in general factor F{H, 9) depends only on the component of the magnetic 
field perpendicular to the interface. 

In conclusion, we have considered the influence of the classical magnetic field on the 
diffusive anomaly in the I-V characteristic of the tunnel junction. We have shown that the 
nonlinear differential tunneling conductance strongly depends on the orientation of the mag- 
netic field with respect to the interface plane. This phenomenon is explained by the different 




(19) 



As it follows from Eq. (|T^, there is a strong anisotropy of the diffusive anomaly with respect 



character of charge spreading after tunnehng for different orientations of the magnetic field. 
The result of this work explains the observation of the recent experimentl. 
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FIGURES 

FIG. 1. The tunnel geometry. The magnetic field vector constitutes angle 9 with the interface 
plane. 

FIG. 2. Diagrams which contribute in the lowest order to two-particle Green's function Eq.(4) 
are shown. Diagrams a) and c) describe correspondingly exchange and Hartree corrections to 
one-particle Green's function in each lead. Diagram b) describes the interaction of electrons in 
different leads. Wavy line is dynamically screened Coulomb interaction. Diffusion vertexes are 
shown by shaded regions. 
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